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Abstract 

We present the general properties of the far-ultraviolet (FUV; 1370-1710 A) continuum background 
over most of the sky, obtained with the Spectroscopy of Plasma Evolution from Astrophysical Radiation 
instrument (SPEAR, also known as FIMS), flown aboard the STSAT-1 satellite mission. We find that 
the diffuse FUV continuum intensity is well correlated with iVm, 100 yum, and Ha intensities but 
anti-correlated with soft X-ray intensity. The correlation of the diffuse background with the direct 
stellar flux is weaker than the correlation with other parameters. The continuum spectra are relatively 
flat. However, a weak softening of the FUV spectra toward some sight lines, mostly at high Galactic 
latitudes, is found not only in direct-stellar but also in diffuse background spectra. The diffuse 
background is relatively softer than the direct stellar spectrum. We also find that the diffuse FUV 
background averaged over the sky has a bit softer spectrum compared to direct stellar radiation. A 
map of the ratio of 1370-1520 A to 1560-1710 A band intensity shows that the sky is divided into 
roughly two parts. However, this map shows a lot of patchy structures on small scales. The spatial 
variation of the hardness ratio seems to be largely determined by the longitudinal distribution of 
OB-type stars in the Galactic plane. A correlation of the hardness ratio with the FUV intensity at 
high intensities is found but an anti-correlation at low intensities. We also find evidence that the FUV 
intensity distribution is log-normal in nature. 

Subject headings: diffuse radiation — dust, extinction — ISM: structure — radiative transfer - 
scattering — ultraviolet: ISM 



1. INTRODUCTION 

Measurement of the interstellar radiation field (ISRF), 
comprised of the direct-stellar and diffuse background 
radiation fields, is of considerable interest because the 
ISRF, in particular at the far-ultraviolet (FUV) wave- 
band, controls the physics and chemistry of the inter- 
stellar medium (ISM). The FUV ISRF is often expressed 
by a simple approximation referred to as the Habing ra- 
diation field (Habing 1968) and was measured for the 
first time by Henry et al. (1977) with two Geiger coun- 
ters on board an Aerobee rocket. Draine (1978) obtained 
another simple analytic expression that appeared to be 
in good agreement with all the previous results at that 
time, and it has been referred to as the Draine's "stan- 
dard" UV background. Gondhalekar et al. (1980) com- 
puted a theoretical model and tabulated the average UV 
ISRF scaled to fit the radiation field observed by the 
S2/68 telescope aboard the TD-1 satellite. Mathis et al. 
(1983) redetermined the ISRF by fitting the spectrum 
derived by Gondhalekar et al. (1980) to power-law and 
scaling up by 15% to bring it into better agreement with 
the measurement of Henry et al. (1980) and extended the 

1 Korea Astronomy and Space Science Institute, Daejeon, Re- 
public of Korea, 305-348; email: kiseon@kasi.re.kr 

2 University of Science and Technology, Daejeon, Republic of 
Korea, 305-350 

3 Space Sciences Laboratory, University of California, Berke- 
ley CA, USA 94720 

Ritter Astrophysical Research Center, University of Toledo, 
Toledo, OH, USA 43606 

5 Korea Advanced Institude of Science and Technology, Dae- 
jeon, Republic of Korea 305-701 

6 Korea Intellectual Property Office, Daejeon, Republic of Ko- 
rea 302-701 



calculation to cover wavelengths from 0.09 ^m to 1000 
/iin. 

The majority of the diffuse FUV background radia- 
tion or the diffuse galactic light is generally believed to 
be of Galactic origin, starlight scattered off by interstel- 
lar dust. Observations of the diffuse background then 
give us information on the scattering properties of the 
dust grains. However, early attempts to measure the dif- 
fuse FUV background were motivated by the hope that a 
large fraction of the radiation in this waveband would be 
extragalactic (see reviews by Paresce & Jakobsen 1980, 
Bowyer 1991, and Henry 1991). There had been seri- 
ous disagreements among observers regarding the spatial 
and spectral variabilities of the FUV background radi- 
ation. Discrepancies among early measurements of the 
diffuse FUV radiation have lead to different interpreta- 
tions as to its dominant contributor, either extragalac- 
tic light or galactic light scattered off dust grains. A 
correlation between the FUV background and Galactic 
neutral hydrogen column density (Nm) was presented 
for the first time by Paresce et al. (1980), thereby im- 
plying that most of the FUV background is Galactic 
in origin, consistent with scattering of the integrated 
FUV radiation field by interstellar dust. Substantial evi- 
dence supporting the conclusion has been reported later 
(Maucherat-Joubert et al. 1980; Joubert et al. 1983; 
Jakobsen et al. 1984; Onaka & Kodaira 1991). 

Some studies, however, have suggested that the cor- 
relation with Galactic latitude is more significant than 
that with Galactic H I column density (Weller 1983; 
Fix et al. 1989), or that the FUV background inten- 
sity shows only a weak correlation with Galactic observ- 
ables (Murthy et al. 1989, 1990). However, Perault et al. 
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(1991) confirmed the correlation with Galactic H I col- 
umn density after subtraction of the inverse sin|6| depen- 
dence from the data obtained with the D2B-Aura satel- 
lite. Hurwitz et al. (1991) found that their data taken 
with the Berkeley UVX spectrometer was three times 
better in dependence with Galactic H I column den- 
sity than with Galactic latitude. Jakobsen et al. (1987), 
Perault et al. (1991), and Sasseen & Deharveng (1996) 
showed also that a correlation exists between the dif- 
fuse infrared background intensity at 100 /zm as mea- 
sured by the Infrared Astronomical Satellite (IRAS) and 
the diffuse background intensity in the FUV wavelengths. 
Sasseen et al. (1995) demonstrated a similarity between 
the spatial power spectra of the IRAS 100 ^m cirrus im- 
ages and the diffuse FUV background images obtained 
with the Far Ultraviolet Space Telescope (FAUST). In 
the scale of a single cloud, Haikala et al. (1995) discov- 
ered a Galactic cirrus cloud, G251.2+73.3, near the north 
Galactic pole in FUV waveband and obtained a good 
correlation between the FUV intensity and IRAS 100 
^m surface brightness. The existence of these proper- 
ties thus provides strong observational support for the 
dust-scattered origin of most, if not all, of the FUV back- 
ground. 

Different contributors to the diffuse FUV background 
have been suggested as well. Duley & Williams (1980) 
proposed fluorescent H2 emission as a significant con- 
tributor to the observed FUV background in addition to 
scattering by interstellar dust grains. Jakobsen (1982) 
showed that the reported correlations could be explained 
by a combination of starlight scattered by interstel- 
lar dust and the H2 fluorescence. Later, the H2 fluo- 
rescence emission lines from the diffuse ISM were ob- 
served with the Berkeley UVX Shuttle Spectrometer 
(Martin et al. 1990) and then by SPEAR (Lee et al. 
2006, 2008; Ryu et al. 2006). Those measurements make 
it clear that the molecular hydrogen emissions are not 
the dominant contributor of the diffuse FUV background. 
Jakobsen & Paresce (1981) suggested that emission lines 
from a hot (> 10 5 K) Galactic coronal gas would con- 
tribute to the diffuse FUV background, although this 
source would not be of sufficient intensity to influence the 
measurements of diffuse FUV continuum background. 
Deharveng et al. (1982) estimated the possible contribu- 
tion of a warm (~ 10 4 K) ionized medium (WIM) to the 
FUV background and attributed a FUV continuum in- 
tensity of ~ 5 — 70 photons cm -2 s _1 sr _1 A -1 (contin- 
uum unit, hereafter CU) to the two-photon continuum 
emission from the WIM (see also Hurwitz et al. 1991). 
Reynolds (1992) showed that the two-photon emission 
from the WIM accounts for approximately 20% of the 
diffuse FUV intensity near 1600 A. 

Disagreements on the spectral shape of the diffuse FUV 
continuum background are also noticeable. Henry et al. 

(1978) measured the diffuse FUV background at the 
north and south Galactic pole regions, using a spectrom- 
eter flown on Apollo 17, and reported a spectrum that 
is flat between 1300 and 1525 A and declining in in- 
tensity toward 1625 A. The observations by Hua et al. 

(1979) and Zvereva et al. (1982) supported the decrease 
of intensity with increasing wavelength at high latitudes, 
using a spectrometer on board the Soviet spacecraft 
"Prognoz-6". They also reported the flat spectrum in 



the intermediate latitudes and a rapid increase of the in- 
tensity with wavelength near and inside the Milky Way. 
However, Anderson et al. (1979) reported some evidence 
for a sharp rise in intensity longward of 1680 A with 
instruments flown on an Aries sounding rocket. The 
data observed by Henry et al. (1980) also showed a rise 
longward of about 1500 A. Murthy et al. (1989) found a 
spectrally flat background over the wavelength band from 
1200 to 1700 A with the Johns Hopkins UVX experiment 
on board the Space Shuttle. The spectrum obtained at 
a high Galactic latitude with the Berkeley UVX FUV 
spectrometer is relatively flat between 1400 and 1850 A 
(Martin & Bowyer 1990). 

Such discrepancies would be better understood with 
spectrophotometric observations over large area of the 
sky. However, only a few studies have mapped a 
substantial fraction of sky in the FUV waveband. 
Gondhalekar et al. (1980) presented the spatial distri- 
bution of the FUV stellar radiation at four passbands 
centered at 1565 A, 1965 A, 2365 A, and 2740 A 
from observations made with the S2/68 sky-survey tele- 
scope. The first coarse map of the diffuse FUV radi- 
ation was presented by Lequeux (1990) and later up- 
dated by Perault et al. (1991) using the data obtained 
with the ELZ spectrophotometer on board the D2B-Aura 
satellite (see also, Maucherat-Joubert et al. 1978, 1980: 
Joubert et al. 1983). The most recent FUV band map 
is from the Narrowband Ultraviolet Imaging Experiment 
for Wide-Field Surveys (NUVIEWS; Schiminovich et al. 
2001), which mapped one-quarter of the sky in FUV radi- 
ation (A = 1740 A, AA = 100 A FWHM) with - 5' - 10' 
imaging resolution, including features such as the Upper 
Scorpius region. 

Spectroscopy of Plasma Evolution from Astrophys- 
ical Radiation instruments (SPEAR, also known as 
Far Ultraviolet Imaging Spectrograph or FIMS) pro- 
vided the first large spectral sky survey of dif- 
fuse FUV radiation (Edelstein et al. 2006a,b). The mis- 
sion observed ~ 80% of the sky and conducted deep 
pointed observations toward selected targets, such as 
supernova remnants (Seon et al. 2006; Nishikida et al. 
2006; Shinn et al. 2006, 2007; Kim et al. 2007, 2010a,b), 
superbubbles (Kregenow et al. 2006; Park et al. 2007; 
Ryu et al. 2008), photoionized H II regions (Park et al. 
2010), and H 2 molecular clouds (Lee et al. 2006, 2008; 
Ryu et al. 2006; Park et al. 2009). The emission lines 
from multiply-ionized carbon, oxygen, silicon, and alu- 
minum toward the north ecliptic and Galactic pole re- 
gions observed with SPEAR/FIMS were presented in de- 
tail by Korpela et al. (2006) and Welsh et al. (2007). 

In this paper, we describe some general results on 
the diffuse FUV (1370-1710 A) continuum background 
observed with SPEAR/FIMS. Emission line sky-maps 
obtained with SPEAR/FIMS is described in Kregenow 
(2007) and Edelstein et al. (2011, in preparation). This 
paper is organized as follows. In §2 we describe the ob- 
servation and data reduction. The FUV intensity maps, 
comparison of the diffuse FUV background with a model 
of the stellar radiation field are shown in §3. Spectral 
shapes of the diffuse FUV background is described in §4. 
We compare the diffuse FUV continuum map with maps 
of 2V H i, Ha, 100 /Ltm, and soft X-ray (1/4 keV) in §5. A 
discussion is presented in §6. Finally, we conclude with 
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a summary of results in §7. 

2. OBSERVATION AND DATA REDUCTION 
2.1. Observation 

SPEAR/FIMS 7 is a dual-channel FUV imaging spec- 
trograph ("Short" wavelength channel [S-band] 900-1150 
A, "Long" wavelength channel [L-band] 1350-1750 A; 
A/AA - 550) with a large field of view (S-band, 4° x 4'. 6; 
L-band, 7°. 4 x 4'.3), designed to observe diffuse FUV 
emission lines. SPEAR is the primary payload on the 
first Korean science and technology satellite, STSAT-1, 
and was launched into a ~ 700 km Sun-synchronous orbit 
on 2003 September 27. The SPEAR/FIMS mission, the 
instruments, their on-orbit performance, and the basic 
processing of the instrument data are described in detail 
by Edelstein et al. (2006a,b). The spectrographs consist 
of the gratings, order baffling, the detector, a shutter- 
slit unit, a mirror unit with field baffling, and a deploy- 
able dust cover. The shutter can be selected to admit 
0%, 1%, 10%, or 100% of the available light for safe and 
photometric observations of bright sources. The instru- 
ment provides for the first time accurate determination 
of the FUV background spectral map covering a sub- 
stantial fraction of sky. The first Galactic map of diffuse 
FUV (L-band) background measured by SPEAR/FIMS 
was presented in Edelstein et al. (2006a). Because S- 
band lacks the Lya-rejection filter that L-band has, the 
major continuum source for S-band is the instrumental 
background (Edelstein et al. 2006b; Kregenow 2007). In 
addition, S-band has lower sensitivity than L-band, as 
well as a higher non-astrophysical background, including 
a strong Ly/3 (1026 A) airglow line. Therefore, in the 
present study, we analyze only the L-band data. Given 
the wavelength dependence of the dust albedo and the 
fact that the composite stellar spectrum, which is dom- 
inated by OB-type stars, peaks near L-band, we expect 
that the astrophysical continuum in S-band is no higher 
than in L-band. 

SPEAR/FIMS observed the sky during the all-sky sur- 
vey observation mode by scanning the entire sky along 
the short axis of the slit, i.e., along the ~ 5' field direction 
of the ~ 8° x 5' FOV. In each orbit, the FOV scanned 
~ 180°, from an ecliptic pole to the opposite ecliptic 
pole, during the eclipse time (~ 25 min), and the pole- 
to-pole scanning drifts 360° along the ecliptic equator for 
one year because of the properties of a sun-synchronous 
orbit. The FUV sky maps for L-band were created by 
binning photon and exposure events using the HEALPix 
tessellation scheme (Gorski et al. 2005) with ~ 1° pixels 
(corresponding to the resolution parameter A s id = 64), 
unless otherwise specified, and L-band wavelength bins 
of 1.0 A. Only the data from 1370 to 1710 A, excluding 
the intense O I airglow line at 1356 A, were used for the 
present analysis. The data were selected for attitude- 
knowledge quality by using times when the derived atti- 
tude error is < 30' and contain 4.1 x 10 7 and 2.0 x 10 7 
photons before and after removal of locally bright stars, 
respectively. The diffuse background map 8 was obtained 
by eliminating photon and exposure events in the locally 
intense pixels (i.e., stars) from the starting sky map with 

7 http://spear.ssl.berkeley.edu/. 

8 The maps presented in the present study can be downloaded 
from the following URL: http://kiseon.kasi.re.kr/FIMS.html. 



~ 14' pixels (A S idc = 256). Each orbital scan was made 
into an image and convolved with an edge detection ker- 
nel shaped like the stellar point spread function. Any 
region with counts 3er above the local background was 
considered a star. All photon events within a 3' radius 
of the star centroid were flagged in the data. This strat- 
egy effectively removed most faint and all bright stars 
(Edelstein et al. 2011, in preparation). The remaining 
photon and exposure events were binned into a map with 
~ 1° pixels (A sidc = 64). 

We repaired orbit number mismatches that were found 
after the publication of Edelstein et al. (2006a,b), and 
included about 100 more orbits other than those used 
in Edelstein et al. (2006a). Sky exposure was derived 
by generating time-marked exposure events and mapping 
the events to a sky position as detailed in Edelstein et al. 
(2006a). The corresponding sky exposure map for the L- 
band with ~ 1° pixels (A s id e = 64), shown in Figure 1, 
covers ~ 80% of the sky and includes features such as 
deep exposures (> 10 ks) toward calibration and pointed 
study fields. About 15% of these observations were taken 
using the 10% shutter position. The exposures taken 
with the 10% shutter aperture were scaled by 0.1, as 
determined in §2.2, and added to the map taken with 
the 100% shutter position. In the figure, the gray color 
represents locations with zero exposure. 

In L-band, there are several spectral features, such as 
atomic emission lines (Si IV A1398, Si II* A1532, C IV 
AA1548, 1551, and Al II A1671), quasi-band structures of 
H2 fluorescence lines, and two-photon continuum emis- 
sion. Since the physical processes that give rise to the 
emission lines are different than dominant continuum 
scattering, the variation of the emission lines may be 
independent from the dust-scattering processes that pro- 
vide the bulk of the signal. If this emission is included 
when estimating the "continuum" intensity, the ability to 
compare the true FUV continuum background and other 
ISM tracers is weakened, as is any discussion of the sta- 
tistical properties of the FUV continuum. To obtain the 
true intensity of the FUV continuum, we may have to 
model all possible emission lines and continuum using 
spectral fits, as in Korpela et al. (2006). However, the 
signal-to-noise ratio of the data with ~ 1° pixels is typi- 
cally only ~ 10 (as described in §2.3), which is too small 
to provide reliable fits. To obtain reliable fits, large pixel 
sizes of ~ 2° — 7.5° (A s ;d c = 32 — 8) were required, espe- 
cially at high Galactic latitudes, depending on locations 
in the sky (see also Edelstein et al. 2011, in prepara- 
tion). In this paper, instead of using spectral fits and 
variable pixel sizes, we excluded the wavelength regions 
of the strongest emission lines (Si II*, C IV, and Al II, 
as shown in Figures 11 and 12) and averaged the data 
in the wavelength regions of 1310-1520 A, 1560-1660 A, 
and 1680-1710 A. Unless otherwise specified, the FUV 
intensities in this paper were obtained in this way. 

We report how the emission lines contribute to the 
total intensity of L-band. Edelstein et al. (2011, in 
preparation) showed that, on average, all atomic emis- 
sion lines and H2 fluorescence lines contribute about 5.1% 
and 2.9% (respectively) to the intensity of SPEAR/FIMS 
L-band. Our independent modeling using spectral fits in- 
dicated that the maximum contribution of emission lines 
is about 10%. Notably, the present results are not signif- 
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icantly different from results estimated using the whole 
L-band, which implies that the emission lines do not 
play a significant role in analyzing the general proper- 
ties of the FUV continuum background. We therefore 
use the terms "diffuse FUV background" and "FUV con- 
tinuum background" interchangeably. The contribution 
of two-photon emission to the diffuse FUV continuum 
background is discussed in §6.3. 

2.2. Instrumental Background 

The SPEAR/FIMS data are affected by instrumental 
(or detector) background due to cosmic rays, radioactive 
decay within the detector, and thermal charged particles 
entering the instrument. These background sources are 
relatively constant in time and uniform across the face 
of the detectors, implying that the instrumental back- 
ground is spectrally flat intrinsically (Korpela et al. 2006; 
Kregenow 2007). 

Since the instrumental background is independent of 
the size of the shutter aperture while the astrophysical 
signal is not, the background can be determined by com- 
paring the data observed with both 100% and 10% shut- 
ter apertures toward the same sky. We thus determine 
the instrumental background rate assuming that the in- 
tensities measured with 100% and 10% shutter aper- 
tures are 7 100% = 7 sky + 7 dot and I w % = a7 sky + 7 dc t = 
a hoo% + (1 - a Kdct, respectively. Here, a, 7 sky , and 7 de t 
denote the scale factor between two shutter apertures, 
the sky intensity measured with 100% shutter aperture, 
and the detector background, respectively. 

We first estimate the detector background intensity av- 
eraged over the whole L-band instead of applying the 
method to every spectral bin, and then calculate the 
instrumental background spectrum in terms of CU. We 
note that each pixel of a HEALpix map is hierarchically 
subdivided into many smaller pixels as the resolution pa- 
rameter A^gide increases. Some of the subdivided pixels 
from a non-zero-exposure pixel at A s ;d c = 64 can have 
zero exposures, especially when the pixel is located at the 
boundary of each orbital scan. Hence, when comparing 
the 100% and 10% maps at A s ;d c = 64, two coinciding 
pixels may not have the same sky coverage at higher res- 
olutions. This discrepancy can produce a large variation 
when comparing 100% and 10% maps with N s id e = 64. 
Thus, we only used pixels (A s ;d e = 64) with 100% sky 
coverage at the level of A s id e = 256, i.e., those pixels 
whose subdivided pixels at the level of N sidc = 256 all 
have non-zero exposures. As shown in Figure 2, over 
the whole L-band, the averaged intensity of instrumental 
background that best fits the measurements is ~ 190± 10 
CU. Similar results that are consistent within ~ 10% 
were obtained using the higher or lower resolution param- 
eters. In the fit, we fully adopted errors in both axes and 
found the best-fit linear straight line between the two in- 
tensities (Press et al. 1992). Because the noise-property 
of the SPEAR/FIMS data is dominated by photon-noise, 
we used only photon-noise in the present study. The esti- 
mated instrumental background was also justified by an 
arbitrary increase or decrease of the background intensity 
and through visual inspection of the combined map. In- 
creasing or decreasing 7d e t by more than ~ 50 CU caused 
the observed sky regions with the 10% aperture to ap- 
pear unrealistically bright or faint in the combined map. 
The scale factor between the two shutter apertures is 



~ 0.10 ± 0.01 as originally designed. 

Denoting the instrumental background rate by i?det 
counts s _1 A -1 in the detector space, the instrumental 
background spectrum Idet(Xi) in units of CU in an i th 
spectral bin is given by 
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where A(Xi), Af2, A A, and AT denote the effective area 
at the wavelength A^, solid angle, wavelength binsize, and 
exposure time, respectively. Applying the maximum like- 
lihood method to Poisson data, the mean instrumental 
background intensity over the L-band is given by 
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where Nt,- m is the number of wavelength bins. We should 
note here that 7d c t ^ Ej 7iet(Ai)/Abi n - Thus the instru- 
mental background spectrum is obtained by 
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and is shown in Figure 3. In the figure, la error limits 
are also shown in dotted lines. The spectrum of the in- 
strumental background is inversely proportional to the 
effective area. 

Using the solid angle of L-band Aft of ~ 1.6 x 
10~ 4 sr and average effective area HjA(Xj)/Nbi n of 
~ 0.19 cm~ 2 , we obtain the detector background rate 
R det = I det E j A(X j ) AO,/ N hi n ~ 0.006 counts s^ 1 A" 1 . 
The instrumental background was independently esti- 
mated to be 0.02-0.04 counts s _1 A _1 by summing many 
shutter-closed dark exposures of 42 ks from 250 orbits 
(Edelstein et al. 2006b). Lee et al. (2006) adopted the 
same method and found a bit weaker detector dark 
background of about 0.01 counts s _1 A _1 . We note 
that the instrumental background rate estimated in the 
present study is even smaller than the value obtained by 
Lee et al. (2006). The discrepancy is attributable to the 
light leakage in shutter-closed dark exposures. 

2.3. Basic Statistics 

The histogram of the sky exposure is shown in Figure 
4(a) on a logarithmic scale. The dotted line in the figure 
shows the exposure distribution for the data observed 
with the 10% shutter aperture only. Figure 4(b) shows 
the intensity distributions of the observed total FUV in- 
tensity (7puv) including both diffuse (7pjjy se ) and direct 
stellar (7 F uy ar ) intensities, and the diffuse FUV back- 
ground only. The FUV background map has median and 
mean intensities of ~ 1250 and of ~ 3170 CU, respec- 
tively. The reason for the difference between the median 
and mean values are obvious from the asymmetry of the 
intensity histogram. The difference in the median value 
from Edelstein et al. (2006a) is mainly due to subtraction 
of the instrumental background and inclusion of more or- 
bits. The total FUV intensity map has median and mean 
intensities of ~ 1610 and <~ 7080 CU, respectively. The 
distribution of signal-to-noise ratio is also shown in Fig- 
ure 4(c). For ~ 1° (A sidc = 64) pixels, about 84% and 
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87% of pixels have significance of 3er or higher from the 
diffuse FUV sky and the total sky, respectively. Mean 
and median signal-to-noise ratios of ipuv 86 are ~ 13 and 
~ 10, respectively, and the ratios of I F °^jy are ~ 17 and 
~ 10, respectively. The statistical properties of the FUV 
intensities estimated over the whole L-band including the 
strong emission line regions are not significantly differ- 
ent. 

3. FUV INTENSITY 
3.1. Morphology 

The maps of the total intensity (ipuv) an£ l the 
diffuse background intensity (/pjfv sc ) observed with 
SPEAR/FIMS are shown in Figure 5 and 6 together 
with the signal-to-noise ratio maps. As observed by 
Henry et al. (1977) and Gondhalekar et al. (1980), the 
sky is strongly anisotropic in FUV. Not only is there 
a strong variation with Galactic latitude, but also with 
longitude. For example, much more FUV radiation is ob- 
served between longitude 180° and 360° than between 0° 
and 180° in the Galactic plane. Gould's belt, inclined at 
~ 18° to the Galactic plane with the direction of tilt to- 
ward Orion, shows up well in Figure 5 (Stothers & Frogel 
1974; Elias et al. 2006 and references therein). 

Obviously, the diffuse FUV background seems to fol- 
low direct starlight in general. However, it will be shown 
in §3.3 and §5 that ipuv sc correlates better with dust, 
traced by 100 /im emission, and H I column density 
rather than with the direct stellar intensity. The largest 
diffuse intensity is toward the Galactic plane and other 
regions where bright early-type stars coexist with suit- 
able columns of interstellar dust, such as the obvious 
feature of the Sco-Cen association and the Vela region. 
We note that the sky is brightest at the Vela region in 
both the total and diffuse background intensities. 

We confirm the existence of a significant depression in 
both the diffuse and total FUV maps at high latitudes 
above less intense segments of the Galactic plane from I = 
20° to I = 60°, as noted in Schiminovich et al. (2001) and 
Gondhalekar et al. (1980). Gondhalekar et al. (1980) 
showed that this direction corresponds to one of the 
most heavily reddened directions in the sky by comparing 
the FUV map with the reddening data of Nandy et al. 
(1978). We also confirmed the strong reddening to- 
ward this direction with the more recent results given 
by Arenouetal. (1992), Chen et al. (1998), and Joshi 
(2005). The dark region in 15° < I < 40° and -6° < b < 
+20° is in fact a molecular cloud referred to as Aquila 
Rift, which emits strong CO line emission and has a max- 
imum extinction of Ay ~ 3 (Dame & Thaddeus 1985; 
Straizys et al. 2003). 

As can be seen in the orthographic projections of Fig- 
ure 6(b), the map reveals lots of detectable structures 
everywhere in the sky. A full examination of the rich de- 
tail in the diffuse FUV map observed by SPEAR/FIMS 
is beyond the scope of this paper. However, it would 
be worthwhile to note that a large structure centered 
at (l,b) ~ (316°, 51°) is a feature produced by dust- 
scattered off starlight from a B1III/IV + B2V spectro- 
scopic binary a Vir with a distance of 87 pc (Park et al. 
2010). An H II region surrounding the star has been 
detected by Reynolds (1985). 



3.2. Dependencies on Galactic Coordinates 

Figures 7 and 8 examine the distribution of /pjfy sc as a 
function of Galactic coordinates. The top panel of each 
figure plots a two-dimensional histogram in the space 
of intensity versus coordinates, while the bottom panel 
shows the median and the standard deviation of intensi- 
ties from the median value within each of coordinate bins. 
The sharp cutoff at low intensities in Figure 4(b) and the 
shape of the latitude profile of the diffuse FUV emission 
in Figure 8 resemble the corresponding features shown 
in the diffuse Ha emission (Figures 10 and 11, respec- 
tively, in Haffner et al. 2003). These similarities between 
the diffuse FUV and Ha emissions may be explained to 
some extent by a simple plane-parallel model. In the ide- 
alized model, the cosecant effect is removed by plotting 
/fuv 6 sin |6| versus sin |6|. This is examined in Figure 9. 
The top panel again plots a two-dimensional histogram 
and the bottom panel shows the median and standard de- 
viation, with the trend for the northern Galactic hemi- 
sphere split from that of the southern hemisphere. It 
should be noted, however, that the plane-parallel model 
turns out to be only approximate in describing the diffuse 
FUV background, as is also the case for the H I gas and 
the WIM (Dickey & Lockman 1990; Haffner et al. 2003). 
The horizontal line in Figure 9 delineates the overall me- 
dian value of I^y sc sin |6| = 525.4 CU. The decrease in 
/diffuse sm 1 6 | at sin | & | < o.05 (i.e., |6| < 3°) is attributed 
to a strong dust extinction in the Galactic disk. 

The diffuse FUV intensity decreases from sin |6| = 0.2 
through 1.0 in both Galactic hemispheres. Similar ef- 
fects have been observed in the Ha and H I distribu- 
tions, and the departure from the csc|& law was at- 
tributed at least in some extent to the presence of the 
Local Bubble (Cox & Reynolds 1987), a localized region 
of very low-density gas around the Sun, as pointed out 
by Dickey & Lockman (1990) for H I. Since the Local 
Bubble is believed to be elongated perpendicular to the 
Galactic plane (Lallement et al. 2003), a larger fraction 
of column density of the ISM is carved out at high lat- 
itudes, yielding the apparent departure from the csc|6| 
law. However, the main reason for the departure in the 
diffuse FUV background is due to the larger contribu- 
tion of dust-scattered FUV radiation at low latitude, as 
shown for the Ha radiation by Wood & Reynolds (1999). 

3.3. Comparison with Stellar Radiation Fields 

Henry (1977) showed that simple addition of the pre- 
dicted UV light from all of the stars listed in common 
star catalogs provides an estimate of the expected ISRF 
in the UV band which is in good agreement with the more 
complex models (see also, Henry et al. 1977). Thus, 
we made a model of direct starlight using the Hippar- 
cos catalog (Perryman et al. 1997; van Leeuwen 2007) 
for stellar location, distance, spectral type, and bright- 
ness. The Hipparcos catalog has a limiting magnitude 
of V ~ 12.4 and was chosen for the calculation of three- 
dimensional locations of stars. The model is basically the 
same as in Henry & Murthy (1993), Henry (2002), and 
Sujatha et al. (2004). 

Based on the spectral type of each star, we derived a 
temperature and effective gravity using the tables from 
Straizys & Kuriliene (1981) and calculated a spectral en- 
ergy distribution for each star with a new grid of Ku- 
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rucz models (Castelli & Kurucz 2003), rather than us- 
ing the original grid of Kurucz models. The intrinsic 
FUV luminosity of each star was obtained by scaling the 
model flux with the factor derived from the formula, as 
in Sujatha et al. (2004): 



AitR 2 
And 2 



F v e- Tv = (jV) 10- y / 2 - 5 , 



(4) 



where R 2 is the scaling factor corresponding to the stellar 
radius, d the distance from the Sun, (fv)o = 3.631 x 10~ 9 

ergs cm" 2 sec -1 A 1 the absolute visual magnitude cor- 
responding to zero magnitude (Bessell et al. 1998), Fy 

the model flux (ergs cnr 2 sec -1 A ) at V band filter, 
V the observed Johnson magnitude, and ry the opti- 
cal depth at V band filter which can be estimated from 
Ty = Ay/1.0863. Here, the absolute visual magnitude is 
calculated from the standard relation Ay = RyE(B-V) 
and Ry = 3.1, and the color excess E(B — V) obtained 
from the observed B — V and the theoretical B — V as 
derived from the Kurucz model spectrum. In many cases 
multiple or complex spectral types were associated with 
a single star; in such cases, we used the first listed spec- 
tral type. We then integrated the stellar flux attenuated 
by the dust distribution along given lines of sight. The 
resulting stellar fluxes in every pixel are summed up to 
yield the "stellar equivalent diffuse intensity (SEDI)" as 
defined by Hurwitz et al. (1991): 



all stars 



rHipparcos \ * 

J SEDI — 



— ^ And 2 f^Healpi: 



(5) 



where Li is the luminosity of each star at the 1370-1710 
A band (excluding the strongest emission line regions) in 
photons s" 1 A" 1 , di is the distance to the star in centime- 
ters, n is the optical depth to each star at the waveband, 
and f^Hcaipix = 47r/./Vpi x is the solid angle (steradians) of 
a pixel. The optical depth at the FUV wavelength is cal- 
culated by rpuv = 2.86(ijfxjv /ay)E(B ~V) where ctfuv 
and ay are the dust extinction cross-sections at the FUV 
wavelength and the V band filter for the average Milky 
Way dust with Ry = 3.1 as from Weingartner & Draine 
(2001) and Draine (2003). We used the resolution param- 
eter -/Vgide = 32 corresponding to an angular resolution of 
~ 1.8° and the number of pixels N p - lx = 12, 288 in order 
to increase the number of pixels including stars. 

We also calculated the SEDI using the TD-1 stellar 
catalog (Thompson et al. 1978). The TD-1 satellite per- 
formed the first UV all-sky survey and measured the ab- 
solute UV fluxes of point sources down to the tenth visual 
magnitude for unreddened early B stars. The TD-1 cat- 
alog presents the stellar fluxes in four bands (each 310- 
330 A wide, centered at 1565, 1965, 2365, and 2740 A) 
for 31,215 stars. The SEDI of the TD-1 stars is defined 

by ^sedi 1 = SiFj/fiHcaipix usin S the stellar fluxes Fj at 
the 1565 A band. Schiminovich et al. (2001) also used 
the TD-1 star catalog to model the diffuse FUV back- 
ground, and they found that stars with 10 times higher 
flux than the TD-1 catalog cutoff contribute about 21 
CU to the FUV background. Therefore, the TD-1 cat- 
alog is complete enough for the analysis of the diffuse 
FUV background. 
First of all, the TD-1 SEDI is compared with the Hip- 



parcos SEDI and the direct stellar intensity obtained 
from SPEAR/FIMS data in Figures 10(a) and (b), re- 
spectively. The Hipparcos SEDI predicts the TD-1 SEDI 
very well in the statistical sense. However, the Hipparcos 
SEDI is a bit higher than the TD-1 SEDI at intensities 
lower than ~ 100 CU, implying that the stars in the TD- 
1 catalog are not complete down to the intensity level 
that was estimated with the Hipparcos catalog, although 
the incompleteness is not so significant. Larger scatters 
in Figure 10(b) than those shown in Figure 10(a) are 
mainly due to the coarse algorithm used to identify stars 
in the SPEAR/FIMS data analysis. However, the cor- 
relation between the TD-1 SEDI and the SPEAR/FIMS 
/|uy ar is good enough to move into the following argu- 
ments. 

We note in Figures 10(c) and (d) that /pw 50 snows 
more or less a general correlation with the Hipparcos 
and TD-1 SEDIs. However, as noted by Hurwitz et al. 
(1991), as long as one avoids the brightest early-type 
stars in the sky by more than a few degrees, there should 
not be a strong correlation between line-of-sight starlight 
and line-of-sight dust illumination. In fact, we found that 
the correlation with the direct stellar intensity is rather 
weaker than that with the infrared 100 /im emission and 
A^hi, as discussed in §5. It is also obvious that the I F uv 
is brighter than both SEDIs, as in Figures 10(e) and 
(f), especially at relatively low intensities I^jy < 10, 000 
CU. The deviation of from the Hipparcos and TD-1 
SEDIs increases as the SEDIs decrease. This trend is a 



diffuse 
FUV 



direct result of the intensity distribution of Jj 
be seen in Figures 10(c) and (d), in which /puy sc shows 
no strong correlation with the SEDIs. More specifically, 
within each "pixel" (~ 1.8°), there is an excess of dif- 
fuse FUV, which is above the estimate from the stellar 
fluxes that are directly detected and/or modeled within 
that pixel and does not correlate with the direct stellar 
flux in the pixel. For low-intensity pixels, the excess of 
diffuse FUV would mainly be due to dust-scattered light 
from stars that are not located along the line of sight 
of the pixels. The conclusion is discussed much more 
clearly in §4.2, which attempts to link a harder spectrum 
at low intensities to this excess. In fact, luminous UV 
stars produce extended dust-scattered halos that are far 
greater than the volume sample of the pixels here (e.g., 
Murhty & Henry 2011). We thus conclude that simply 
adding stellar flux underestimates ISRF, especially at 
the low intensity of Jisrf ;$ 10, 000 CU, against the re- 
sults of Henry (1977) and Henry et al. (1977). Three- 
dimensional radiative transfer models may help to as- 
sert what portion of the ISRF is the result of the dust- 
scattered FUV background. 

It would be worthwhile to note that the contribution 
of unresolved faint stars to our diffuse FUV background 
map is negligible. The diffuse background or the total 
intensity would be correlated at low intensity pixels with 
the model SEDI if this was the case. However, it is ob- 
vious in Figure 10(c) and (d) that this is not the case. 
The difference of more than an order between Isedi and 
the mean value of ipuv a ^ ^SEDi ~ 50 CU in Figures 
10(e) and (f) is unlikely due to the contribution from un- 
resolved stars. We thus conclude that residual starlight 
is not a significant contributor to the diffuse FUV back- 
ground. The contribution from faint stars seen in Fig- 
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ure 10(a) is much smaller than the difference between 
-?sedi and /puv m Figures 10(e) and (f). Many studies 
(i.e., Hurwitz et al. 1991; Henry 2002) have estimated 
the contribution of unresolved faint stars with several 
methods in their analysis of the FUV background data, 
and they all have concluded that the contribution is not 
significant. Additional support for this conclusion will 
be provided in §4 from an argument based on hardness 
ratio. We also note that the large difference between the 
diffuse FUV intensity and the model SEDIs cannot be 
due to H2 fluorescence lines and two-photon continuum 
emission, which may account for ~ 2.9% (Edelstein et 
al. 2011, in preparation) and ~ 4 — 9% (see §6.3) of the 
diffuse intensity, respectively. 

Zodiacal light, sunlight scattered by interplanetary 
dust, is also negligible at these wavelengths since the so- 
lar GO spectrum falls steeply at the FUV wavelengths 
(Tennyson et al. 1988), and we have confirmed that 
our data show no systematic trend with ecliptic lat- 
itude. Scattering of diffuse Lya radiation (geocoro- 
nal/interplanetary) as determined from laboratory mea- 
surements is also negligible. The insignificance of un- 
resolved faint stars, zodiacal light, and Lya scatter- 
ing is also justified from the agreement of the so-called 
"isotropic component" estimated in §5 with the previous 
measurements . 

4. FUV SPECTRUM 
4.1. Spectral Variations 

The direct stellar and diffuse background intensities 
should show strong anisotropy not only in intensity but 
also in spectral shape. In Figures 11 and 12, the av- 
erage spectra of the diffuse background, direct starlight, 
and total radiation in various Galactic longitude and lat- 
itude ranges, respectively, are shown. The direct stellar 
spectra were obtained by subtracting the diffuse back- 
ground spectra from the total spectra in each coordinate 
range. We note that the FUV intensities are strongest in 
240° < I < 270° and in -10° < b < 0°, where the Vela 
region is located, as was seen in Figure 6. The Vela re- 
gion is the brightest not only in emission lines due to the 
Vela supernova remnant but also in the FUV continuum 
(cf. Nishikida et al. 2006). 

The average background intensity is about the same as 
the direct stellar intensity in all coordinate ranges shown 
in the figures, and the total intensity is approximately 
two times higher than the direct stellar intensity. Ignor- 
ing the latitude range b < —10°, where exposures and 
signal-to-noise ratios are small, the intensity ratio of the 
diffuse background to the direct starlight decreases with 
Galactic latitude, and then increases at the highest lati- 
tude range of 70° < b < 90°. The trend is easily under- 
standable, at least qualitatively, by adopting the results 
from a simplified spherical dust model with an embedded 
source. The scattered to absorbed stellar flux ratio in the 
model increases monotonically with optical depth and 
can easily be of the order of unity or higher for optical 
depth of r = 1 - 2 (Witt et al. 1982). The optical depth 
of t - 1.26 at 1550 A is estimated at N m = 1 x 10 21 
cm -2 using the standard gas-to-dust ratio (Bohlin et al. 
1978) and the Milky Way dust extinction cross-section 
(Weingartner & Draine 2001; Draine 2003). Therefore, 
the approximate equality between the diffuse and direct 



stellar intensities at the Galactic plane is obtained. The 
decrease of the diffuse to direct stellar intensity ratio with 
Galactic latitude is also explained since the optical depth 
decreases with latitude. However, at the highest Galactic 
latitude range 70° < b < 90°, the "isotropic component" 
of ~ 300 CU discussed in §6 contributes significantly to 
the diffuse background in addition to the scattered radi- 
ation of Galactic starlight in the latitude range, and thus 
causes the ratio to increase. 

A weak rise in the diffuse spectra longward of about 
1550 A is shown in most of the latitudes in Figure 12, 
although the FUV continuum spectra are flat in gen- 
eral. However, the dependence of the rise on the Galac- 
tic latitude is not clear. The rise is also noticeable at 
120° < I < 180° in Figure 11, where most of observa- 
tions were performed at high latitudes as seen in Figure 
1. It should be noted as well that the softer the diffuse 
spectrum is, the softer the direct stellar spectrum also 
is. A much sharper rise longward of ~ 1550 A than that 
found here was obtained by Henry et al. (1980), in which 
they concluded that the rise was due to very large num- 
bers of faint stars contributing to the ISRF. Later, Henry 
(2002) reanalyzed the data and attributed the rise to 
an actual rise in the interstellar grain albedo. However, 
these features may not be due to a single cause. Firstly, 
a stronger extinction at a shorter wavelength leads to the 
reddening or softening not only of direct stellar, but also 
of diffuse spectra, unless the dust albedo decreases with 
wavelength as will be shown in §6 (Figure 25). Secondly, 
if there is no significant contribution from bright early- 
type stars near a given sight line, faint late-type stars 
should contribute, at least in part, to the spectra in ad- 
dition to the dust-scattered light from distant early-type 
stars. 

It should also, more importantly, be noted that the dif- 
fuse background spectrum is always softer than the direct 
stellar spectrum. The rise of the dust albedo with wave- 
length is in fact closely related to the relative softness 
of the scattered to the direct stellar spectra, as will be 
discussed in §6. We also note that both the direct stellar 
and diffuse spectra are harder, or bluer, in the Galactic 
longitude range 240° < I < 270°, where the Vela re- 
gion dominates the radiation field, than other longitude 
ranges, implying the existence of lots of hot stars and/or 
less extinction in the region. A relatively sharp rise long- 
ward of about 1690 A at |6| > 50° is also found mostly 
in the diffuse background spectrum. However, the origin 
of this rise is not clear at this time. 

We now compare the spatially-averaged total (di- 
rect stellar + diffuse) FUV spectrum observed with 
SPEAR/FIMS with the previous models and observa- 
tions in Figure 13. It can be easily noted that all pre- 
vious results are higher, up to a factor of more than 3 
for the Draine's model, than the SPEAR/FIMS data. In 
Figure 12, the observed total spectrum is stronger only 
at -10° < b < 0° than the Draine's "standard" radia- 
tion field, while the intensities in other latitude ranges 
are in general weaker. The relative weakness of the ob- 
served data compared to the model ISRFs indicates that 
the previous works are mostly biased toward bright re- 
gions in Galactic plane. In Figure 13, we also plotted 
the average spectra of the diffuse background and the 
direct starlight. The figure shows once again that the 
direct starlight and the diffuse background contributes 
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approximately the same amounts to the total spectrum 
on average, implying underestimation of the ISRF when 
only direct starlights are added. 

The spectral shapes of the model ISRFs calculated 
by Gondhalekar et al. (1980) and the Habing field are 
more or less similar to the SPEAR/FIMS data, while 
the Draine's "standard" and Mathis background radia- 
tion fields are a bit harder than that of SPEAR/FIMS. 
The observed spectrum of Henry et al. (1980) is softer 
than the SPEAR/FIMS data. We also note that the 
diffuse background spectrum averaged over all the sky 
is a bit softer than the stellar spectrum, and that dis- 
crete spectral features such as atomic emission lines are 
mainly shown in the diffuse spectrum only, as already 
shown in Figures 11 and 12. The relative softness or 
redness of the diffuse background compared to the direct 
stellar spectrum is once again attributed to rise of the 
dust-scattering albedo at a longer wavelength, as will be 
discussed in §6. 

4.2. Hardness Ratio 

In order to further investigate spectral variations of 
the diffuse background, we plotted a map of the hard- 
ness ratio defined as the ratio of average intensities at 
AA1370 - 1520 A and at AA1560 - 1710 A (excluding 
AA1660 - 1680 A to avoid the Al II emission line) in 
Figure 14. In the figure, we selected a color scheme to 
represent the region with relatively blue or hard spec- 
tra as blue colored. The spectral hardness map shows 
highly patched structures and no clear trend along with 
Galactic latitudes. Instead, we note that the sky can 
be largely divided into two blocks mainly along Galactic 
longitudes. In fact, we smoothed the hardness ratio map 
using a gaussian kernel with full width at half maximum 
of ~ 3.5° and saturated the color range to clarify the divi- 
sion of the sky based on the hardness ratio. A bit harder 
radiation is observed between longitude 180° and 360°, 
where stronger radiation is found in Figure 6, than be- 
tween 0° and 180°. Interestingly, we found that the hard- 
ness ratio of the diffuse background follows in general 
the spatial distribution of OB-type stars in the Galac- 
tic plane. In Figure 15, we plotted maps for numbers of 
OB-type stars and of A- type stars using the TD-1 stellar 
catalog. We note from the figure that the Galactic plane 
can be divided into three regions: 240° < I < 360° where 
OB-type stars are predominantly found, 30° < I < 90° 
and 180° < I < 240° where A-type stars are dominant, 
and 90° < I < 150° where both OB- and A- types are 
rare. In the Galactic longitudes where OB-type stars are 
dominant, the hardness ratio is generally higher than the 
other regions. 

We also plotted the ratio versus Galactic longitude and 
latitude in Figures 16 and 17, respectively. In Figure 
16, the median value of the hardness ratio at each lon- 
gitude bin is generally lower in 0° < I < 180° than in 
180° < / < 360°. On the other hand, Figure 17 shows 
no clear trend along Galactic latitude. In order to con- 
firm the close relation between the hardness ratio and 
the number of OB-type stars, we calculated the average 
hardness ratio and the number of OB-type stars within 
each of the 3° longitude bins, and plotted the average 
hardness ratio versus the number of OB-type stars in 
Figure 18. The figure shows a clear correlation between 
the average hardness ratio and the number of OB-type 



stars. We also note a dip at 120° < I < 180° in Figure 
16. The dip does not seem to be an artifact caused by 
a relatively low exposure at low latitudes in this longi- 
tude range. In the light of the argument described in the 
previous paragraph, the softness in the longitude region 
would be a direct consequence of no bright stars in the 
range. 

Figure 19 shows the hardness ratio versus the diffuse 
FUV intensity. Rises of the hardness ratio at high and 
low intensities are noted in the figure. To understand 
the trend, we calculated the hardness ratio versus the 
direct stellar intensity using the previous stellar model; 
the result is shown in Figure 20. In the calculation of 
the hardness ratio, we incorporated the wavelength de- 
pendence of dust extinction as well. We note that the 
model ratios predict the observed values nicely at high 
intensities (i|{f{T > 2000 CU), of which most of the 
radiation is from low Galactic latitudes. Therefore, we 
attribute the rise of the hardness ratio at high intensities 
(/diffuse > 2000 CU) mostly to the spectral hardness of 
stars that are located at low Galactic latitudes. 

However, the slight rise of the hardness ratio at low 
intensities (/puy 86 < 2000 CU) cannot be explained 
by stellar types at high Galactic latitudes. One possi- 
ble explanation for this rise is the contribution of dust- 
scattered starlight from earlier stars rather distant than 
nearby stars that are located along the same line of sight 
for each pixel. This is supported by Figures 10(c) and 
(d) , which show clear excesses of the diffuse background 
above the nearby direct stellar intensity. In fact, our 
preliminary results of the Monte Carlo simulation, which 
will be presented in detail elsewhere, show that B-type 
stars are significant sources of the diffuse FUV back- 
ground even at the high Galactic latitudes where most 
stars are A-types. Henry (2002) also noted that B-type 
stars contribute the largest amount to the ISRF in his 
simplified model. We thus conclude that both the rise 
of the hardness ratio and the excess of the diffuse back- 
ground over the direct stellar intensity at low intensities 
are largely due to the contribution of B-type starlight 
scattered off dust grains. However, at the lowest inten- 
sities of less than a few hundred CU, the "isotropic com- 
ponent" would contribute significantly. A detailed study 
on the relative significance between the two components 
may need extensive radiative transfer models that are 
beyond the scope of this paper. 

5. COMPARISON WITH VARIOUS GALACTIC 
QUANTITIES 

We now examine correlations of the observed FUV 
background with various galactic quantities, not only 
with the H I column density and 100 fim emission but 
also with the Ha and soft X-ray (1/4 keV) emissions. 
Although extensive correlation studies between the dif- 
fuse FUV background, the H I column density, and 100 
/im emission have been performed, no detailed correla- 
tion study with other waveband observations has been 
done. The all-sky map of H I column density are from the 
Leiden/Argentine/Bonn (LAB) Survey (Kalberla et al. 
2005). The 100 fim and Ha emission maps are from 
Schlegel et al. (1998) and Finkbeiner (2003), respec- 
tively. The soft X-ray (1/4 keV) map is from the ROSAT 
All-Sky Survey (RASS) map (Snowden et al. 1997). The 
H I data interpolated onto a HEALPix projection were 
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obtained from the Legacy Archive for Microwave Back- 
ground Data Analysis (LAMBDA). 

Correlations between the diffuse FUV intensity and 
other quantities are shown in Figure 21 in log-log scales. 
Two-dimensional histogram images and contours are 
shown in the figure instead of the usual pixel-by-pixel 
comparison to clarify their correlation. Correlation coef- 
ficients estimated in logarithmic scale are also shown in 
Figures 21. It is also worthwhile to note that the cor- 
relation coefficient between the diffuse FUV background 
and the direct stellar radiation is ~ 0.46, much less than 
those for TVhi and Iioo/im- In Figure 22. we examine 
the correlations in linear-linear scales, especially at high 
Galactic latitudes (b > 25°). The general similarity in 
the latitude dependence of the data sets (the H I column 
density, 100 /im, Ha, and FUV emissions) may dominate 
their correlations. To reduce the latitude dependence, we 
therefore calculated the average intensities multiplied by 
sin |6| within each of the A sin |6| = 0.02 latitude inter- 
vals. The correlations between the quantities multiplied 
by sin \b\ are shown in Figure 23, in which only the data 
sets observed at high Galactic latitudes (|6| > 30°) were 
examined, in order to minimize the dust extinction effect 
at low Galactic latitudes. The numbers outside and inside 
the parentheses in Figure 23 are correlation coefficients 
calculated in linear and logarithmic scales, respectively. 
Here, we note that Figure 23 shows strong correlations 
between emission tracers, implying that the correlations 
between the FUV background and other emission trac- 
ers are not simply due to the latitude dependence of the 
Galactic observables. 

5.1. Correlation with iVm and lOOfim emission 

Correlation plots of the diffuse FUV background with 
H I column density and 100 /im intensity are shown 
in Figures 21(a) and (b), respectively. There is large 
scatter in the correlation plots. We note that the 
diffuse FUV continuum intensity saturates at about 
Ahi ~ 10 21 cm -2 , which was found for the first time 
by Hurwitz et al. (1991). The same effect of saturation is 
also shown in comparison with the 100 /im intensity. This 
saturation at a high intensity can be attributed to strong 
dust absorption of the FUV radiation at low Galactic 
latitudes. Flattening at a low intensity, which has been 
often attributed to extragalactic background, is also seen 
in Figures 21(a) and (b). 

As is obvious in the figures, the diffuse FUV intensity 
for a given Nui or 100 /im brightness varies more than by 
a factor of 2-3, up to an order of magnitude, with direc- 
tions in the sky. Witt et al. (1997) argued that most of 
the observed variation is a result of the spatially, severely 
anisotropic, ISRF. Since the diffuse FUV background de- 
pends not only on the interstellar dust but also on the in- 
situ stellar radiation field, the usual pixel-by-pixel com- 
parison might be affected by the strong anisotropic stellar 
radiation field. We thus averaged the FUV data corre- 
sponding to a given H I column density or a given 100 
/im brightness and compared the average FUV intensity 
with TVhi and 100 /im brightness. In this way, the effect 
of the spatially varying stellar radiation field may be av- 
eraged, and only the effect due to H I or dust column 
density can be investigated. The average FUV intensity 
versus TVhi and Iioo^m are shown in Figure 22, together 
with the best regression lines. In the figure, the data 



used in the fit are denoted by filled circles, and the data 
not used in the fit by hollow circles. In the fit, we used 
only the data for b > 25°. The best-fit regression lines 
are given by 

Ifuv c = (1-49 ± 0.07) x 1f 5 HI „ + (271.2 ± 26.7) CU, 

(6) 

and 

^fuv c = (158.3±11.7)x /l00 ^ m +(243.1±44.4) CU (7) 
MJy/sr 

for H I column density and 100 /im emission, respectively. 

Paresce et al. (1980) reported variations in the slopes 
and brightness axis intercepts corresponding to the var- 
ious scans. As noted by Jakobsen et al. (1984), there 
is no unique, canonical relationship between background 
intensity and a hydrogen column density that is obeyed 
everywhere in the sky. Principally, this result is not sur- 
prising, since the detailed viewing geometry and the illu- 
minating UV stellar radiation field are highly anisotropic, 
and several of the physical parameters, such as dust 
properties, are expected to vary from region to region 
in the ISM. However, it is obvious, as shown in Figure 
21(a) and (b), that there is a general correlation in the 
global scale between the diffuse FUV background and 
both the neutral hydrogen column density Nm and IR 
100 /im emission. The linear correlation has been mea- 
sured by many different investigators, and the results 
published before 1991 were summarized in Bowyer (1991) 
as I^f c ~ (0.3-2.5) x (A m /10 18 cm- 2 )+300 CU, which 
is consistent with our result. 

The flaring and large scatters shown in these linear- 
linear correlation plots are in general due to large varia- 
tions of the diffuse FUV intensity as noted in the log-log 
correlation plots of Figure 21. Note that the larger scat- 
ters are found at the higher mean intensities, the prop- 
erty that will be discussed in more detail in §6. We also 
note a straight-line correlation for Jioo^m < 8 MJy/sr 
in Figure 22(b). Witt et al. (2008) found that there 
was also a linear correlation between the optical sur- 
face brightness of high-latitude cirrus clouds and their 
100 /im intensity at least up to ~ 8 MJy/sr. This is 
the range in which these objects remain optically thin 
and single scattering by dust grains is dominant. The 
systematic deviation from a straight-line correlation for 
hoOfim > 8 MJy/sr in Figure 22(b) would be due to 
the fact that at the higher Jioo/xm values we are entering 
lower galactic latitudes where the effects of the strongly 
forward-directed scattering phase function of dust pro- 
duces strongly enhanced FUV intensities compared to 
those produced by the scattering at mostly larger scat- 
tering angles at higher latitudes, where the 100 /im in- 
tensities are generally low. 

5.2. Correlation with the Diffuse Ha Emission 

In linear scale, the best regression line between i|j? v so 
and /net is given by 

/diffuse = ( 456 ± 23.5) x I^L + (309.2 ± 31.4) CU, (8) 

R 

where R denotes Rayleigh (1 R = 10 6 /47r photons cm -2 
s^ 1 sr _1 ). The average /puv se versus 1h and best-fit 
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line are shown in Figure 22(d). Only the data for b > 25° 
were used in the fit. 

We note a very nice, unexpected, correlation between 
the FUV background and the diffuse Ha emission, some 
of which is known to be emitted from the WIM. As can 
be seen in Figures 21 and 22, the low intensity behavior 
in the correlation of -fpuy 30 with Jh„ is slightly different 
from in the correlation with either A'hi or Iioo/xm- Figure 
21 shows that the correlation between /p}yy sc and Jn tt is 
overall straight down to the lowest intensities, while other 
relations become flat or curved a bit upward. Moreover, 
in Figure 22, at the low intensity of I^ sc (~ 400 - 500 
CU), the diffuse FUV intensity seems to decouple from 
the 100 /jm intensity (and a bit from the H I column 
density), while the correlation between /p^fy 36 and In a 
is still retained. These results seem to indicate that the 
FUV background has a stronger correlation with the Ha 
emission, at least qualitatively, than with the other emis- 
sion tracers, although the correlation coefficient between 
7pjjy sc and In a is not significantly higher than the other 
correlation coefficients (Figure 21). 

Figure 23 shows the results after the removal of 
inverse-sin |6| dependence due to the plane-parallel 
medium, offering clearer evidence that the correlation 
between /puy 50 and /hq is stronger than other cor- 
relations. Since the far-IR emission from interstellar 
dust correlates well with neutral hydrogen emission (e.g., 
Boulanger & Perault 1988), the diffuse FUV continuum 
correlates well with both iioo^m and Nm at the same sig- 
nificance, as indicated by the correlation coefficients in 
Figure 23(a) and (b). If the correlation between the FUV 
and Ha backgrounds is simply due to the co-existence 
of dust and Ha-emitting gas, the correlation of /puy sc 
and Iho would have the same significance as with either 
Jioo/jm or Nui. On the other hand, the correlation with 
the Ha background would be enhanced when the physical 
origin of a large portion of the Ha background is analo- 
gous to the FUV background. We may then conjecture 
that the better correlation of the FUV with Ha intensi- 
ties is attributed to the similarity of the physical origins 
of the diffuse FUV and Ha emissions. The fact that the 
correlation between /pjfv sc and Irq appears to exist not 
only on a large scale but also in smaller-cloud scales sup- 
ports the idea wherein their physical origins are at least 
in part related. We discuss this connection briefly in §6.1 
and in more detail in a paper in preparation. 

We here consider several possible explanations for the 
relationship between the FUV continuum and Ha emis- 
sions, assuming that the physical origins of the back- 
grounds are indeed closely related with each other. First 
of all, a large fraction of the correlation between /puy sc 
and In a may be caused by the two-photon continuum 
emission in FUV wavelengths, which originates from the 
WIM. However, the contribution of two-photon emission 
to the FUV continuum {I^/Ina ~ 4 — 9%, as found in 
§6.3) is not large enough to explain the correlation rela- 
tion Eq. (8). 

Wood & Reynolds (1999) estimated that the scattered 
Ha intensity at high Galactic latitudes is ~ 5 — 20% of the 
total Ha intensity. Recently, Mattila et al. (2007) and 
Lehtinen et al. (2010) found that in some dust clouds, 
the Ha radiation is mostly due to scattering by dust 
grains that are illuminated by a Ha-emitting source off 



the line of sight. Witt et al. (2010) showed that a sub- 
stantial fraction of the diffuse high-latitude Ha back- 
ground is caused by interstellar dust scattering of Ha 
photons that originate elsewhere in the Galaxy. In par- 
ticular, they found many extended regions where the 
fraction of the scattered Ha intensity is of an order of 
50% or higher. Witt et al. (2010) and Dong & Draine 
(2011) also found that the most likely average fraction of 
scattered Ha intensity is about 20% of the total Ha inten- 
sity. Thus, the scattered component of the Ha emission 
may strengthen the correlation between the diffuse Ha 
and FUV intensities, especially at high latitudes. 

It should be noted here that the differences between 
the data in Figure 22(b) and Figure 22(d) accord with 
the increasing significance of the dust-scattered Ha emis- 
sion with latitude. At low latitudes the diffuse Ha 
emission is mostly from ionized gas and not related to 
dust-scattering, while the FUV background is from dust- 
scattering. Therefore, at low latitudes and high values 
of the diffuse FUV intensity, a much wider range of vari- 
ability in the Ha intensity is found for a given value of 
the diffuse FUV intensity (Figure 22(d)) while the dif- 
fuse FUV intensity still correlates with 100 /im emis- 
sion showing an enhancement due to the effects of strong 
forward-scattering (see §5.1). On the other hand, at high 
latitudes and low values of the diffuse FUV intensity, 
much of the structure in the faint Ha background may 
be due to scattering. Thus, we are seeing a linear cor- 
relation with the diffuse FUV intensity there. In other 
words, the Ha contribution from dust-scattering becomes 
increasingly important at higher latitudes. Such a strong 
linear correlation seems to imply that a significant por- 
tion of Ha at high latitudes is due to dust-scattering. 
A more detail comparison of the diffuse FUV and Ha 
maps in both large and small scales is needed to address 
how large portion of the diffuse Ha emission is due to 
dust-scattering. 

Finally, there is one possibility beyond the scattered 
Ha that is worth considering. The ISM may have low- 
density paths and voids that allow for ionizing photons 
from midplane OB stars to reach and ionize gas many 
kiloparsecs above the Galactic plane (Wood et al. 2010). 
Pathways that provide lower than average densities to 
high latitudes could enhance both scattered FUV and 
gas ionized by the Lyman continuum from the same OB 
stars that are much closer to the Galactic plane. Com- 
bined with the scattered Ha emission, the presence of 
such pathways could provide a stronger correlated link 
between FUV and Ha. In this regard, we note that if 
ionizing photons could use low-density paths to pene- 
trate a turbulent ISM and reach high latitude clouds, 
such as LDN 1780, the resulting ionization would be lim- 
ited to a thin outer shell and the ambient medium of the 
clouds (Witt et al. 2010). Wood et al. (2010) estimated 
that only -0.05 pc thick skin of LDN 1780 would be 
photoionized. For such morphology, the optically thick 
core regions would appear relatively dark, which contra- 
dicts the observations described in Mattila et al. (2007) 
and Witt et al. (2010). However, we note that this does 
not preclude ionizing photons travelling to high latitudes. 
Photoionization would contribute to a portion of the Ha 
intensity seen in LDN 1780, although dust-scattering is 
the dominant process for the Ha intensity on the cloud. 
Therefore, further investigation is needed to disentangle 
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which process (dust-scattered Ha emission versus in-situ 
ionized gas) is more significant in explaining the correla- 
tion between the diffuse FUV and Ha emissions. 

A weak saturation effect in the correlation between 
7puy so and ijja > which is similar to but much weaker than 

those in the correlation of /fuv s ° vers us Am and hoo/im, 
is also noticeable. The weak saturation seen in the cor- 
relation with /hq may be easily understood through the 
fact that the WIM can also be approximated with a 
plane-parallel model but the extinction cross-section at 
FUV is higher than that of Ha. 

5.3. Correlation with Galactic Latitude 

If the diffuse FUV background correlates primarily 
with Galactic latitude the FUV background inten- 
sity is described by C\ esc |6| + C%. In addition to these, 
we also noted in the previous section that the diffuse 
FUV intensity can be approximated with a simple plane- 
parallel model defined by Cicsc|6|. We thus fitted the 
data with two trial functions, one with two free parame- 
ters and the other with only one free parameter. 

The dependence on esc | & was also obtained by av- 
eraging the diffuse FUV intensity for a given Galactic 
latitude. Best regression functions are then given by 



rdiffusc 
'FUV 



846.7 ±96.1 



and 



sin \b\ 



rdiffusc 
i FUV 



(-457.2 ± 100.5) CU, (9) 



412.3 ± 10.3 
sin |6| 



CU 



(10) 



for two-parameters and one-parameter models, respec- 
tively. The average intensity of the FUV background 
versus \b\ are also shown in Figure 22. The regression 
functions with two and one parameters are denoted by 
solid and dashed lines, respectively. The value of the two- 
parameters function at sin \b\ = 1 is a bit lower than the 
value obtained by Fix et al. (1989), but still within their 
error range. We note that the two-parameters model is 
better in explaining the dependence of /pjfv 80 on Galactic 
latitude. However, the physical interpretation of the two- 
parameters model is problematic, as discussed in §6.4. 

5.4. Comparison with the Soft X-ray 

Comparison of the FUV background with the soft X- 
ray emission has been mentioned by Joubert et al. (1983) 
and Zvereva et al. (1982). Especially Zvereva et al. 

(1982) claimed there was a good correlation of the FUV 
intensity in high latitude range > 30°) with soft X- 
ray brightness. They interpreted the correlation as an 
indication of the significant contribution of the hot ISM 
to the diffuse FUV background. However, Joubert et al. 

(1983) argued that hot gas has too small an emission 
measure to contribute significantly to the FUV back- 
ground. Soft X-ray emission has been observed ex- 
tensively with ROSAT and is well-known to be anti- 
correlated with H I column density. Indeed, we found 
a strong anti-correlation between the diffuse FUV inten- 
sity and the soft X-ray background, as shown in Figure 
21(d). This anti-correlation between the diffuse FUV 
background and the soft X-ray also supports the dust- 
scattered origin of the diffuse FUV background. 



6. DISCUSSION 

6.1. Log-normal nature of the FUV intensity 
distribution? 

In presenting the results, we used logarithmic scales 
to visualize large orders of magnitude in observed FUV 
intensity. Here, we show that the distribution of the dif- 
fuse FUV intensity seems to be well represented by a 
log-normal distribution or, equivalently, a Gaussian dis- 
tribution in a logarithm scale, although its detailed study 
is beyond the scope of the present paper. 

Hill et al. (2008) found that the histogram of the Ha 
emission line observed from the WIM fits well with 
a log-normal distribution. They attributed the log- 
normal nature to the turbulence in the WIM. A log- 
normal distribution is the expected probability density 
function of density and/or column density for the ISM 
with a density structure established by turbulence (e.g., 
Vazquez-Semadeni & Passot 1999). Since the scattered 
FUV intensity is directly related to the dust column den- 
sity, the log-normal property of the distribution of the 
diffuse FUV intensity is a natural consequence of the den- 
sity structure and the turbulence property of the ISM. 
Thus in Figure 24, we plotted intensity histograms of 
7puy so sin \ b\ for various Galactic latitude ranges together 
with best-fit log-normal functions, and found that the 
intensity distributions are indeed well represented with 
log-normal functions. Here, we removed the Magellanic 
Clouds and 19 large-scale regions (the same regions as in 
Hill et al. 2008) with significantly enhanced Ha emission 
to avoid discrete ionized regions where the contribution 
of two-photon emission to /puv sc may be significant. We 
also removed some high intensity "streaks" (most visi- 
ble in the left panel of Figure 6(b)) that show up in the 
exposure maps. 

Log-normal distribution shows an extended tail to 
large intensity values in a linear scale. Some evidence 
of the extended tail in an intensity histogram has al- 
ready been noticed in the previous studies. Joubert et al. 
(1983) noted that the distributions of the points in corre- 
lations between UV brightnesses and H I column density 
are not symmetrical with respect to the linear regression 
lines. There was a high flux tail in the distribution of 
the UV intensities at a given value of A(H I). The inten- 
sity histograms in Figure 2 of Paresce et al. (1980) also 
show such extended tails. These regions with UV inten- 
sity excesses were rejected for the correlation studies in 
Joubert et al. (1983) and Perault et al. (1991) by elimi- 
nating the points in the tail with counts larger than the 
median by several times the dispersion and iterating the 
procedure. Joubert et al. (1983) attributed this prop- 
erty to an excess FUV radiation in certain regions of the 
sky, perhaps due to two photon emission by ionized gas. 
However, Deharveng et al. (1982) and Reynolds (1992) 
showed that the emission from the WIM is unlikely to 
contribute significantly to the general diffuse UV back- 
ground. Perault et al. (1991) found that the excess in 
the intensity histogram does not obey the Poisson distri- 
bution; they attributed the high-count skewness and tail 
to the contribution of the stars brighter than m(UV) = 8 
and unresolved faint stars. However, we argued that con- 
tribution of residual stars to the diffuse FUV background 
is not significant. In addition, the faint late- type stars 
cannot explain the hardness ratio observed at low inten- 
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sity regions, as shown in §4. Schiminovich et al. (2001) 
also noted a large variance in the FUV intensity and at- 
tributed the observed variance to the "cosmic variance" 
predicted by the multi-cloud model combined with multi- 
ple scattering and three-dimensional radiation field vari- 
ance. We thus believe that the anomalous high inten- 
sity tails shown in previous studies are in fact due to 
an intrinsic property of the log-normal nature of the 
diffuse FUV background. A detailed statistical prop- 
erty of the dust density distribution can be investigated 
through a radiative transfer modeling of the starlight 
dust-scattered off in turbulent media. 

It may also be possible to isolate high-intensity tails 
in the histograms and identify extended local features 
in the FUV maps. Some locally enhanced features were 
found, but bright regions generally had a patchy appear- 
ance. We also found that some of the enhanced features 
in the FUV sky coincided with the Ha sky, but not al- 
ways. We note here that the observed intensity toward 
a cloud could be higher, lower, or equal to its surround- 
ings (as demonstrated in Mattila et al. 2007), depend- 
ing on the scattered intensity, incident intensity from be- 
hind the cloud, and cloud optical depth. Therefore, the 
high-intensity tails in histograms do not necessarily cor- 
respond to local dust clouds. In addition, because of the 
difference in dust-extinction cross-sections at the FUV 
wavelengths and Ha line, there would not be a strict 
correlation between diffuse FUV and Ha emissions, even 
when all the FUV and Ha background photons origi- 
nate from dust-scattering. The origin of both the high- 
intensity tails in the histograms and the diffuse Ha line 
emission may be elucidated by more detailed analyses of 
the dust clouds observed in both the FUV wavelengths 
and Ha line, along with self-consistent radiative transfer 
models. But, such analyses are beyond the scope of the 
present paper. 

6.2. Relative softness of the diffuse radiation 

We found a weak reddening of the diffuse radiation rel- 
ative to the direct starlight in §4. In order to understand 
the trend, we calculated the scattered to direct intensity 
ratio for a simplified case where a point source is embed- 
ded in a homogeneous spherical dust cloud. As noted in 
Witt et al. (1982), the analytical expression derived by 
Code (1973) approximates the scattered to the unatten- 
uated luminosity ratio very well. In the approximation, 
the scattered intensity is given by 

-^scatt _ 2 _ T 

~lT ~ (1 + e« T + (1 - e-« T ~ 6 
where ( = y/ (1 — a) / (1 — ag) and £ = 
y(l — a) (1 — ag), and the attenuated direct starlight 
by ^dircc = Loe~ T . Here, the quantities Lq, a, g, and 
t refer to the unattenuated stellar luminosity, grain 
albedo, phase function asymmetry factor, and radial 
optical depth of the spherical cloud measured from the 
center to the outer edge. 

We calculated ratios L^- llec /L 0y L scatt /L a , and 
■^scatt/^dircc as a function of wavelength using the ex- 
tinction cross-section c(A), a(A), and g(X) from Draine 
(2003) for an illustrative case of t(1550 A) = 2 in Fig- 
ure 25. Here, Lq, Ldircc, and L sca tt are the input, direct, 
and scattered luminosities. Note that the albedo a(X) 



is an increasing function with the wavelength. For com- 
parison, we also plotted the results calculated for three 
hypothetical cases: (1) the case of a constant albedo of 
a(1550 A) and a variable asymmetry factor g(X), (2) 
the case where both a and g are constants obtained 
at 1550 A, and (3) the case in which the wavelength- 
dependence of a(X) is reversed so that the albedo de- 
creases with the wavelength. The figure shows that as- 
suming the asymmetry factor to be a constant does not 
significantly affect the result while the albedo does. As is 
evident, the directly-escaped spectrum is redder or softer 
than the input spectrum (Figure 25a). When a constant 
albedo is assumed, we obtain a scattered spectrum red- 
dened compared to the unattenuated input spectrum, 
but bluer than the (attenuated) directly-escaped spec- 
trum because of the rapid increase of the absorption with 
the wavelength. However, as the albedo increases with 
the wavelength in the model of Weingartner & Draine 
(2001) and Draine (2003), the scattered spectra rise much 
more rapidly with the wavelength than the attenuation, 
thereby giving a relatively-red scattered spectrum com- 
pared to the directly escaped spectrum. In the case in 
which the wavelength-dependence of a(A) is reversed, the 
scattered spectrum become bluer relative to both the 
unattenuated input and the reddened directly-escaped 
spectra, which is contrary to the other cases. The same 
trend is observed when the optical depth is varied. We 
therefore conclude that the relative softness of the diffuse 
background radiation compared to the stellar radiation 
field is attributed to the increase of dust albedo with the 
wavelength. 

In §4, we noticed a weak rise in the diffuse spectra 
longward of 1550 A. The rise itself tells nothing about 
the wavelength-dependence of dust albedo unless the in- 
cident unattenuated spectrum Io(X) is known. Instead, 
the relative softness or hardness of the scattered radia- 
tion compared to the directly observed radiation can give 
direct information about the wavelength-dependence of 
dust albedo. In this regard, we note that the dust albedo 
provided by Bruzual et al. (1988) is constant or decreases 
with the wavelength in FUV, and thus is not consistent 
with our results. A more detailed analysis on the spec- 
tral shape can be obtained by radiative transfer models 
including wavelength dependence in calculating stellar 
radiation fields rather than that performed for a single 
wavelength band as have been done by Witt et al. (1997) 
and Schiminovich et al. (2001). 

6.3. Contribution of two-photon continuum emission 

The WIM, which is expected to have a temperature of 
~ 8000 K, must be a source of two-photon continuum in 
the FUV wavelengths, potentially impacting the inter- 
pretation of the FUV continuum background (Reynolds 
1992). The two-photon continuum intensity at ~ 1600 
A due to the WIM with a temperature of 8000 K is es- 
timated to be 7 27 = 57.4(/ H q/R) CU (Reynolds 1990, 
1992; see also Figure 26). Using the average Ha in- 
tensity distribution /jja ~ 1.2 esc |fo| R, which was ob- 
tained from rather limited observations, Reynolds also 
obtained I2-* ~ 70 esc |6| CU. However, using the result 
of Hill et al. (2008) obtained from the full WIM obser- 
vations, the distribution of Ha for \b\ > 10° is fitted on 
average by Jhq ~ 0.625 esc \b\ R (see also Dong & Draine 
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2011). This leads to a lower two-photon intensity of 
I 2l ~ 36 csc \ b\ CU, which accounts for - 9% of the FUV 
background (Eq. 10). 

The two-photon continuum's contribution to the dif- 
fuse FUV intensity is reduced even further when 
the model of WIM developed by Dong & Draine 
(2011) is adopted. Recent observations using 
the Wilkinson Microwave Anisotropy Probe (WMAP) 
have found that the ratio of the free-free ra- 
dio continuum to Ha is surprisingly low in the 
WIM (Davies et al. 2006; Dobler & Finkbeiner 2008; 
Dobler, Draine, & Finkbeiner 2009). Dong & Draine 
(2011) considered a three-component model consisting 
of a mix of (1) photoionized hot gas, (2) gas that is re- 
combining and cooling after removal of a photoionizing 
source, and (3) cool H I gas. In their standard model for 
explaining the observed intensity ratios of the free-free 
radio continuum to Ha and [N II] A6583 to Ha, the au- 
thors assumed that the scattered fraction of the Ha orig- 
inating from the hot gas is /<™ fl) =0.2, while the frac- 
tions of the hot, cooling, and H I gases are = 0.22, 
^(cooling) = Q and ^Hi) = Q respectively The 

temperature of the hot gas was found to be T( hot > = 9100 
K. The cooling gas cooled down rapidly from its initial 
temperature T (cooling ) - 10,000 K to - 400 K (Fig. 3 of 
Dong & Draine 2011). 

Using the two-photon spectral profile (Kwok 2007) , the 
hydrogen recombination coefficient a2S for case B (Table 
3 of Martin 1998), and the Ha emission rate (Eq. 10 
of Dong & Draine 2011), we estimated the two-photon 
intensity averaged over 1370-1710 A from hot and cool- 



ing ionized gas (Figure 26): 7. 



and /. 



(cooling) 
27 



(hot) 
27 



59.1(WR) CU, 



26.0(/ H a/R) CU. The two-photon in- 



tensity of the cooling gas is mostly determined by the 
gas with a temperature of 400 K, because the cooling 
gas is 400 K or lower for more than 98% of its lifetime 
(Fig. 3 of Dong & Draine 2011). The total two-photon 
intensity originating from the ionized gas is then given 

, T _ f (cooling) r (cooling) f (hot) r (hot) _ da[t /-r>\ 
by i2 7 - / H a 1 2 1 + JHq 1 2 1 - 28.4(J H a/R) 

CU. The distribution of two-photon emission is then 
l2 7 ~ 18 csc |& CU, corresponding to ~ 4% of the FUV 
background intensity. 

In summary, two-photon emission may account for ~ 
4 - 9% of the diffuse FUV background. Although the 
two-photon emission may be significant in discrete bright 
H II regions, the contribution of two-photon emission to 
the diffuse FUV background is much smaller than has 
previously been suggested. The large variation of /puy sc 
(from a factor of 2-3 up to an order of magnitude) is also 
unlikely to be the result of two-photon emission. 

6.4. Isotropic component 

The phrase "isotropic component" has been used to re- 
fer to the diffuse intensity extrapolated to Am = 0, or, 
alternatively, to \b\ = 90°, and is frequently assumed 
to be extragalactic background. Wright (1992) claimed 
that the isotropic extragalactic component of the diffuse 
FUV background should be obtained by extrapolating 
to csc(|6|) = 0, instead of to csc |6| = 1 as practiced 
by Fix et al. (1989). Witt & Petersohn (1994) discussed 
systematic errors arising in the method of extrapolat- 



ing high latitude measurements of the diffuse FUV back- 
ground to Ahi = or to csc( |6| ) = to derive the level 
of the cosmic extragalactic background. They argue that 
the derivation of the cosmic background from the extrap- 
olation of the diffuse FUV background to Ahi = is more 
likely to yield a reliable result. They found, in fact, that 
extrapolations of the predicted FUV background inten- 
sities through a Monte-Carlo radiative transfer model to 
csc |6| = produce negative values when it should pro- 
duce a zero intercept. We also found from the SPEAR 
data that the extrapolation to csc |6| = indeed produces 
a negative value of ~ —400 CU. 

A discussion of the origin of the isotropic component is 
beyond the scope of the present paper. Onaka & Kodaira 
(1991) observed intensities as low as 300 CU at high 
Galactic latitudes. Schiminovich et al. (2001) placed 
constraints on the extragalactic background 200 ±100 
CU. We also found the isotropic component consistent 
with previous authors, but a bit higher. It should be 
noted that a similar level of the isotropic component was 
found from the relationship between the diffuse FUV and 
Ha emissions (Eq. 8). However, it is not clear at this 
time that all of the ~ 300 CU, the extrapolation value 
of the average FUV-Ahi correlation to Ahi = 0, cor- 
responds to the extragalactic background. In fact, we 
observed the FUV intensities lower than 300 CU as is ob- 
vious in Figure 4. Witt et al. (1997) analyzed the FUV 
data observed in the FAUST experiment (Bowyer et al. 
1993) and estimated a contribution of about 700±200 
CU that is uncorrelated with Galactic parameters. They 
argued that the extragalactic component of more than 
300 CU is excluded by the lower limit directly measured 
by Onaka & Kodaira (1991) and concluded that most 
likely value of isotropic extragalactic background radi- 
ation would be 160±50 CU, after correction for Galactic 
extinction. 

7. SUMMARY 

In this paper, we reported some general results on the 
diffuse FUV background and the total FUV radiation 
observed with SPEAR/FIMS. Our main conclusions are 
as follows: 

1. Strong anisotropies are found in not only the diffuse 
background intensity but also in the hardness ratio 
(1370-1520 A to 1560-1710 A). 

2. We found a good correlation between the diffuse 
FUV background and other Galactic quantities, 
such as H I column density, 100 /im emission, and 
Ha emission. The correlation of the FUV con- 
tinuum background with the diffuse Ha emission 
seems to be better than with the other quantities. 
Anti-correlation with the soft X-ray was also found. 
Correlation of the FUV background with the direct 
stellar intensity is rather weaker than the correla- 
tion with other waveband observations, such as 100 
/im and Ha emissions. 

3. The fact that a linear correlation of the FUV back- 
ground with the diffuse Ha emission is found at low 
intensities seems to indicate that the dust-scattered 
component of Ha becomes increasingly important 
at higher latitudes. 
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4. The spatially-averaged total FUV intensity, in- 
cluding the direct stellar and diffuse FUV inten- 
sities, observed with SPEAR/FIMS is weaker than 
the well-known ISRF models, such as the Habing, 
Mathis, and Draine models. 

5. The spectrum of the diffuse FUV background is 
in general flat and a bit softer than the directly- 
escaped stellar spectrum. The relative softening 
can be attributed to the rise of dust albedo as the 
wavelength increases. 

6. The hardness ratio seems to follow in general the 
longitudinal distribution of OB-type stars in the 
Galactic plane. 

7. The "isotropic component" obtained by comparing 
the FUV background with other ISM tracers, A/hi, 
100 /im, and Ha emissions is consistent with the 
previous results. 

8. Evidence that the intensity histogram of the diffuse 
background is well represented by log-normal dis- 
tribution is found, as for the diffuse Ha emission. 

The above results well accord with the fact that most 
of the diffuse FUV background is scattered light of the 



FUV stellar radiation by the interstellar dust. Witt et al. 
(1997) and Schiminovich et al. (2001) developed radia- 
tive transfer models to interpret the diffuse FUV back- 
ground. Our future work will deal with a Monte-Carlo 
simulation model incorporating three-dimensional stel- 
lar distribution obtained from the Hipparcos catalog 
and three-dimensional dust distributions modeled by 
several authors such as Drimmel & Spergel (2001) and 
Drimmel et al. (2003). The model compared with the 
SPEAR/FIMS data should provide a better understand- 
ing to the observational results. 
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Fig. 1. — (a) Mollweide and (b) orthographic projections of sky exposure for the SPEAR/FIMS L-band observations. The map was 
made using the HEALPix scheme with ~ 1° pixels (N s a e = 64). The intensity scales are logarithmic across the color bars. Exposure time 
observed with both 100% and 10% shutter apertures are combined, (a) Galactic coordinates centered at ((,6) = (0°,0°) with the longitude 
increasing toward the left are shown with latitude and longitude lines on a 30° grid, (b) The left and right sides of orthographic projections 
are centered at the northern and southern Galactic polar caps, respectively, and their longitude increases clockwise and counter-clockwise, 
respectively. In both projections, I = 0° is at the six o'clock position. 
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Fig. 3. — Instrumental background spectrum. The spectrum is essentially the inverse of the effective area curves, as shown in 
Edelstein et al. (2006b). An la error range for the average of detector background value is also shown with dotted lines. 
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Signal-to-Noise Ratio 

Fig. 4. — Histograms of (a) exposure time (s), (b) the FUV intensity (photons s — 1 cm" 1 sr _1 A^ 1 , CU) over the observed sky, and (c) 
signal-to-noise ratio with bins scaled logarithmically. The histograms are made with about 1° pixels (iV s ide = 64). A background of 190 
CU was removed from the data to take into account the instrumental dark background. 
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(b) Log 10 (I lotal CU) 




Fig. 5. — (a) Mollweide and (b) orthographic projections of the combined total (diffuse + direct stellar) FUV map observed with 
SPEAR/FIMS, after removal of the instrumental background. The left and right sides of orthographic projections are centered at the 
northern and southern Galactic polar caps, respectively, and their longitude increases clockwise and counter-clockwise, respectively. In 
both projections, I = 0° is at the six o'clock position. The intensity scales are logarithmic across the color bars. 





(d) Log 10 ([S/N] total /3) 




-1.0 ^^^^^^^^^^^h 3.0 



Fig. 5. — (c) mollweide and (d) orthographic projections of signal-to-noise ratio of the total FUV map. The left and right sides 
of orthographic projections are centered at the northern and southern Galactic polar caps, respectively, and their longitude increases 
clockwise and counter-clockwise, respectively. In both projections, Z = 0° is at the six o'clock position. The intensity scales are logarithmic 
across the color bars. 
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Fig. 6. — (a) Mollweide and (b) orthographic projections of the diffuse FUV background map observed with SPEAR/FIMS, after removal 
of locally intense pixels (stars) and the instrumental background. The left and right sides of orthographic projections are centered at the 
northern and southern Galactic polar caps, respectively, and their longitude increases clockwise and counter-clockwise, respectively. In 
both projections, I = 0° is at the six o'clock position. The intensity scales are logarithmic across the color bars. 
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(c) Log 10 ([S/N] dlffuse /3) 
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Fig. 6. — (c) mollweide and (d) orthographic projections of signal-to-noise ratio of the diffuse FUV map. The left and right sides of 
orthographic projections are centered at the northern and southern Galactic polar caps, respectively, and their longitude increases clockwise 
and counter-clockwise, respectively. In both projections, I = 0° is at the six o'clock position. The intensity scales are logarithmic across 
the color bars. 
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Fig. 7. — The diffuse FUV intensity versus Galactic longitude. The top panel plots the two-dimensional histogram of the diffuse FUV 
intensity as a function of Galactic longitude. The bottom panel plots the median value within each of the 3° longitude strips from the 
data. The vertical extent of each plotted bar is determined from the average deviation about the median within that bin. 
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Fig. 8. — The diffuse FUV intensity versus Galactic latitude. The top panel plots the two-dimensional histogram of the diffuse FUV 
intensity as a function of Galactic latitude. The bottom panel plots the median value within each of the 2° latitude strips from the data. 
The vertical extent of each plotted bar is determined from the average deviation about the median within that bin. 
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Fig. 9. — ^fuv S ° s ™ H versus sm H- The top panel plots the two-dimensional histogram as a function of Galactic longitude. The bottom 
panel plots the median value within each of the A sin |6| = 0.025 latitude strips from the data. The vertical extent of each plotted bar is 
determined from the average deviation about the median within that bin. The black and red points represent the data observed toward 
northern and southern Galactic latitude, respectively. The horizontal line denotes the median Tpijv sm H va l ue - 
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Fig. 10. — Comparison between the TD-1 stellar equivalent diffuse intensity (SEDI), the Hipparcos SEDI, and the SPEAR/FIMS 
intensities, (a) The TD-1 SEDI versus the Hipparcos SEDI, (b) the direct stellar intensity observed with SPEAR/FIMS versus the TD-1 
SEDI, (c) the diffuse FUV background versus the Hipparcos SEDI, (d) the diffuse FUV background versus the TD-1 SEDI, (e) the total 
(direct stellar + diffuse) FUV intensity versus the Hipparcos SEDI, and (f) the total FUV intensity versus the TD-1 SEDI. Here, the 
resolution parameter N B j^ e = 32 corresponding to angular resolution of ~ 1.8° is used. The TD-1 SEDI was estimated using the 1565 A 
band, and the others were calculated over the same wavelength bands 1370-1520, 1560-1660, and 1680-1720 A. 



OBSERVATION OF THE FUV CONTINUUM BACKGROUND 



27 




1 , , , , , ; 

1400 1450 1500 1550 1600 1650 1700 
Wavelength (A) 



20 


- 240' < ] < 270° 


SI II" C IV 


A] 11 


15 








10 


*" ~ "** **" "* — *" 






5 


















1400 1450 1500 1550 1600 1650 1700 
Wavelength (A) 




1 - i 
01 , , , , , : J 

1400 1450 1500 1550 1600 1650 1700 
Wavelength (A) 



Fig. 11. — Average FUV spectra in various Galactic longitude ranj 
and dashed curves represent the total (direct stellar + diffuse) and t 
also indicated. 
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Fig. 12. — Average FUV spectra in various Galactic latitude ranges. Solid curve shows the diffuse FUV background spectrum. Dotted 
and dashed curves represent the total (direct stellar + diffuse) and the direct stellar spectra, respectively. Some important ionic lines are 
also indicated. 
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Fig. 13. — SPEAR/FIMS spectra averaged over all the observations are compared with the three models of the ISRF. The four models 
are from Habing (1968), Gondhalekar et al. (1980), Draine (1978), and Mathis et al. (1983). The spectrum measured by Henry et al. (1980) 
is also compared. 



30 SEON ET AL. 

( a ) 'l370-1520 / 'l560-1710 




HAD ^ 2J1 



Fig. 14. — (a) Mollweide and (b) orthographic projections of the hardness ratio (1370-1520 A to 1560-1710 A) map. The Al II emission 
line was excluded in estimating the hardness ratio. The left and right sides of orthographic projections are centered at the northern and 
southern Galactic polar caps, respectively, and their longitude increases clockwise and counter-clockwise, respectively. In both projections, 
I = 0° is at the six o'clock position. The scales are linear across the color bars. 
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(b) Number of A— type stars 




Fig. 15. — (a) Number of OB-type stars and (b) number of A-type stars from the TD-1 stellar catalog. The maps were made using the 
resolution parameter -/V s id e = 16 (with ~ 3.7° pixels). 
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Fig. 16. — Hardness ratio (1370—1520 A to 1560—1710 A) versus Galactic Longitude. The Al II emission line was excluded in estimating 
the hardness ratio. 
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Fig. 18. — Correlation between the average hardness ratio (1370-1520 A to 1560-1710 A) and the number of OB-type stars. The average 
hardness ratio and the number of OB-type stars were calculated within each of the 3° longitude strips, and then the average hardness ratio 
versus the number of OB-type stars was plotted. The Al II emission line was excluded in estimating the hardness ratio. 
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Fig. 20. — Stellar model prediction of hardness ratio (1370-1520 A to 1560-1710 A) versus intensity (1370-1520, 1560-1660, and 
1680-1710 A). The Al II emission line was excluded in estimating the hardness ratio. 
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Fig. 21. — Correlation of the diffuse FUV background with (a) neutral hydrogen column density, (b) 100 (im emission, (c) the diffuse Ha 
emission, and (d) soft X-ray (1/4 keV band) background. The left panels show two-dimensional histogram and the right panels contours of 
the histograms. Numbers in the left panels are the correlation coefficients estimated in logarithmic scale. The contours correspond to 0.7, 
0.3, and 0.1 of the maximum values of the histograms. 
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Fig. 22. — The FUV intensity versus (a) neutral hydrogen column density, (b) 100 fim emission, (c) Galactic latitude, and (d) the Ho 
intensity. Best fits in (a), (b), and (d) are lines with (a) slope 1.49 CU/(10 18 cm" 2 ), offset 271 CU, (b) slope 158 CU/(MJy/sr), offset 243 
CU, and (d) slope 456 CU/R, offset 309 CU. Best fit curves in (c) are 847csc|b| - 457 CU for solid curve, and 412csc|6| CU for dashed 
curve. Filled and hollow circles denote the data points, respectively, used and not used for the fit. 
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Fig. 23. — Correlation of the diffuse FUV continuum intensity with (a) neutral hydrogen column density, (b) 100 jim emission, and (c) 
the diffuse Ha intensity after the removal of inverse-sin |6| dependence due to the plane-parallel medium. Data points were obtained from 
latitude ranges of sin|6| > 0.5 (|6| > 30°). Numbers outside and inside the parentheses are correlation coefficients estimated in linear and 
logarithmic scales, respectively. 
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Histograms of 7{|g&J se sin |6| in logarithmic scale. Log-normal fits are also shown. The standard deviations and the mean 



Fig. 24. — 

values of the distribution in logarithmic scale are also shown 
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Fig. 25. — Scattering and absorption in a homogeneous dust medium, (a) The directly-escaped and (b) scattered spectra relative to 
the input luminosity of the point source, (c) The scattered spectrum relative to the attenuated direct stellar spectrum. Here, various 
combinations of the wavelength-dependences of the dust albedo and asymmetry factor are examined. 
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Fig. 26. — Two-photon continuum intensity that accompany 1 R of Ha emission from ionized interstellar hydrogen. The continuum 
intensity is shown for three temperatures: 9100 K, solid line; 8000 K, dashed line; 400 K, dotted line. 



